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NRL/ ARPANET  REMOTE  DATA  COLLECTION, 
PROCESSING,  AND  DISPLAY  TERMINAL 


I . INTRODUCTION 

NRL  is  involved  in  a research  study  to  determine  the  ability  to 
cross  correlate  acoustic  signals  emanating  from  sources  in  the  ocean 
after  propagating  over  long  and  widely  separated  distances. 

During  the  fall  of  1973  over  800  hours  of  digitized  recordings 
were  recorded  at  several  acoustic  stations.  The  ability  to  reduce  and 
analyze  this  data  presents  an  enormous  task  when  the  number  of  computa- 
tions and  the  number  of  parameters  that  can  be  applied  to  each  computation 
are  considered.  In  pursuit  of  solutions  to  these  problems,  extensive 
use  has  been  made  of  the  Illiac  IV  computer  complex  for  high  speed 
computation.  A computer  terminal  system  is  located  at  NRL  for  visual 
and  statistical  interpretation  of  the  Illiac  IV  data  output. 

It  is  the  purpose  of  this  report  to  define  the  relationship  of  the 
Illiac  IV  complex,  the  ARPANET,  NRL' s dedicated  computer  installation 
and  the  hardware/software  considerations  used  to  unify  all  of  these  into 
a productive  system.  Included  in  a separate  report  are  flowcharts  of 
the  system  and  technical  details  of  the  subprograms. 

The  aspects  of  the  total  system  can  be  categorized  into  four  broad 
general  areas. 

A.  ARPA  system  modules 

B.  NRL  system  modules 

C.  NRL-ARPANET  interrelational  modules 

D.  Unifying  software 

Each  of  these  areas  will  be  discussed,  outlining  the  individual 
components  and  'esponsibilities  as  they  relate  to  the  overall  system  for 
accessing  raw  field  data  and  displaying  computed  results. 

System  Overview 

Successful  reduction  of  recorded  data  to  a form  desirable  by  NRL  for 
examination  of  computed  results  requires  many  steps  and  subfunctions 
before  meaningful  parameters  can  be  extracted  from  the  data  base. 

During  the  summer  and  fall  of  1973,  acoustic  signals  were  recorded  at 
the  ARPA  Acoustic  Research  Center  for  later  analysis  by  the  Illiac  IV 
complex. 

Since  all  signals  from  each  remote  site  were  multiplexed  onto  a 
single  magnetic  tape,  they  must  be  separated  and  reformatted  into  a form 
more  suitable  to  the  Illiac  IV.  In  addition,  the  source  signal  that 
was  recorded  on  board  the  ship  must  be  included  into  the  data  base 
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recorded  by  the  Acoustic  Research  Center  personnel.  The  approximate 
starting  time  of  the  acoustic  signal  must  be  located  for  each  site  to 
avoid  lengthy  computations  on  the  Illiac  that  do  not  include  the  signal 
of  interest.  Furthermore,  the  tapes  must  be  adjusted  relative  to  each 
other  since  the  acoustic  signal  arrives  at  each  remote  site  at  different 
times  due  to  the  long  propagation  path  involved.  Finally  these  signals 
must  be  digitally  filtered  and  converted  to  a floating  point  notation 
suitable  to  the  Illiac  instruction  code  and  scaled  according  to  relative 
signal  level  of  the  recorded  data. 

Once  the  tapes  are  prepared  in  the  above  manner,  NRL  personnel  can 
request  that  selected  data  be  transferred  to  the  Illiac  IV  complex  for 
Inclusion  into  the  I4-Tenex  for  analysis  by  the  correlation  algorithm  in 
the  Illiac  IV  processor.  This  correlation  program  was  written  by  the 
staff  of  SCI  and  is  documented  in  references  1,  2,  and  3. 

Fig.  1 is  a system  overview  of  the  various  modules  that  are  included 
in  reduction  of  new  data.  The  heavy  outlined  areas  and  blocks  indicate 
hardware  facilities  of  the  ARPANET,  NRL  and  the  Illiac  IV  computer  complex. 
Dotted  lines  and  blocks  indicate  logical  steps  and  Interactions  performed 
to  completely  reduce  magnetic  tape  input  to  displays  at  NRL  for  analysis 
and  parameter  extraction. 

As  shown  in  step  IB  of  Fig.  1,  a request  for  new  experimental  data  is 
initiated  at  NRL.  This  Information  is  sent  via  the  ARPANET  to  personnel 
at  the  Institute  for  Advanced  Computation  (lAC).  The  experimental  data  is 
read  into  the  disk  memory  and  held  there  temporarily  until  the  correla- 
tion program  can  be  executed  on  the  Illiac  IV  processor.  Upon  proper 
completion  of  the  Illiac  program,  the  input  data  is  generally  deleted  in 
order  to  allow  space  for  new  input  data.  If  it  is  desired  to  vary 
parameters  on  the  same  section  of  data,  the  Tenex  disk  areas  may  be 
preserved  until  completion  of  all  computations. 

If  it  is  desired  to  change  parameters  on  the  correlation  program,  the 
procedure  outlined  in  lA  of  Fig.  1 is  followed.  All  parameters  that  can 
be  adjusted  on  the  correlation  program  are  contained  in  a special  Tenex 
disk  file.  Fig.  2 contains  a complete  listing  of  this  file  and  displays 
all  the  parameters  that  may  be  adjusted  on  the  correlation  program.  This 
parameter  program  will  be  explained  in  more  detail  in  the  following 
sections. 

As  shown  in  lA  of  Fig.l,  the  parameter  file  is  modified  at  NRL.  Before 
these  parameters  can  be  utilized  by  the  Illiac  IV,  however,  they  must  be 
compiled  by  the  Borroughs  B6700  computer;  this  is  accomplished  during 
the  normal  batch  processing  of  a job  submission  to  the  I4-Tenex. 

The  I4-Tenex  is  the  master  control  computer  for  all  job  submissions 
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and  for  control  of  all  other  computers  and  subassemblies.  Requests  for 
new  data  tapes,  modification  of  parameters,  execution  of  the  Illiac 
array  program,  storage  of  resulting  data,  and  transmission  of  results  via 
the  ARPANET  are  all  initiated  or  controlled  by  the  I4-Tenex. 

All  of  these  functions  of  the  I4-Tenex  are  not  automatic.  Each  step 
must  be  controlled  by  constructing  a special  file  of  commands  (called 
ACL)  that  direct  each  process  or  subassembly  in  the  Illiac  IV  complex. 

These  ACL  statements  are  contained  in  a Tenex  disk  file  known  as  the 
Primary  Input  File  (PIF);  they  can  be  modified  via  the  NRL  operator  to 
define  new  input  data,  obtain  new  parameters  from  the  B6700,  assign  disk 
areas  for  computational  results,  etc.  Logically  the  PIF  is  shown  in 
step  2 of  Fig.  1.  The  PIF  is  described  more  fully  in  another  section 
and  is  shown  in  Fig.  3. 

One  of  the  responsibilities  of  the  PIF  is  to  assign  areas  on  the 
Illiac  IV  disk  memory  (I4DM).  This  disk,  logical  unit  3 of  Fig.  1,  is 
associated  with  a particular  NRL  job  only  during  actual  execution  of  the 
correlation  program.  It  should  not  be  confused  with  the  Tenex  disk, 
logical  unit  4 in  Fig.  1,  which  is  used  for  long  term  storage  of  all 
programs  and  input/output  data.  The  I4DM  is  specially  designed  to 
support  the  fast  computation  speed  of  the  parallel  processor.  All  input 
data  (from  magnetic  tape)  must  be  placed  on  the  I4DM  before  the  parallel 
processor  begins.  All  output  computations  are  placed  on  the  I4DM  as  the 
correlation  program  executes  on  the  processor. 

Computational  output,  logical  step  3B  of  Fig.  1,  on  the  I4dm  must  be  j 

transferred  to  the  Tenex  disk  (or  other  bulk  storage),  under  direction  j 

of  the  PIF,  if  it  is  to  be  saved  and  transmitted  to  NRL  for  analysis.  ' 

When  the  PIF  has  completed  the  data  transfer  to  Tenex  disk,  the  I4DM  is 
deallocated  from  the  present  job  and  the  next  user  in  the  batch  queue 

takes  command  of  the  I4DM.  I 

Once  the  output  computations  are  available  on  the  Tenex  disk  memory, 
the  data  may  be  reformatted  for  graphical  display,  statistical  printout, 
etc.,  and  transmitted  to  NRL.  The  programs  for  transmission  reside  in 
the  l4-Tenex  and  the  Tenex  disk  data  base  is  accessed  according  to 
directions  supplied  by  the  NRL  operator.  These  interactive  commands 
allow  the  operator  to  select,  via  the  PDP-11  terminal  system,  the  mode  of 
transmission  and  the  output  device  at  NRL  to  receive  the  data.  Examples 
of  presently  available  outputs  are  shown  in  Figs.  5B,  6,  7,  8,  9 and 
explained  in  following  sections.  The  PDP-11  terminal  system  is  shewn 
in  Fig.  10. 

Major  emphasis  is  being  placed  on  preparing  a more  complete  operating 
system  for  the  PDP-11  system.  This  system  will  permit  automated  inter- 
action, under  operator  control,  with  the  data  base  and  allow  the  operator 
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to  flexibly  examine  the  computed  results  and  parameters  by  various 
methods. 


Introduction  of  New  Correlation  Parameters 


The  actual  submission  of  an  executable  job  requires  several  pre- 
liminary steps  before  actual  execution  on  the  parallel  processor  can 
begin.  The  first  operation  is  the  evaluation  of  preliminary  data  at  NRL 
to  determine  which  parameters  are  to  be  varied  for  a given  test  run.  If 
the  data  had  been  run  before,  adjustments  to  the  parameters  may  be 
required  to  provide  further  studies  and  analyses  of  the  signals. 

The  number  of  options  available  in  the  program  is  quite  large. 
However,  under  normal  conditions  only  a small  number  will  be  changed  at 
any  one  execution  on  the  Illiac  processor.  Under  a typical  run,  those 
parameters  needing  adjustment  would  be  modified  by  a program  called 
PARMVO. SOURCE. 

A typical  listing  of  this  program  is  shown  in  Fig.  2.  In  line  1,  a 
frequency  of  33  Hz  is  defined  by  the  parameter  F.  (Note  that  <:=>  is 
the  symbol  for  definition  and  <;>  is  a delimiter  that  separates  one 
parameter  from  another.) 

LOWE STLEVFMT  is  a program  control  parameter  that  limits  the  formatted 
(suitable  for  typing)  output.  Depending  on  the  application  only  higher 
level  integration  times  may  be  needed  for  the  present  analysis.  (The 
meaning  of  these  integration  times  is  defined  more  fully  in  another 
section).  This  parameter  is  defined  as  an  integer  whose  value  is  associa- 
ted with  the  integration  time.  Typically  this  value  is  set  to  5 so  that 
only  the  integration  levels  of  5 and  above  are  output  from  the  program 
and  sent  to  permanent  storage.  Seven  is  the  highest  level  integration 
time  permitted  by  the  present  program.  If  this  parameter  is  defined  as 
8 only  summary  statistics  are  output.  A value  greater  than  8 suppresses 
all  output.  It  should  be  emphasized  that  these  outputs  are  suitable 
only  for  outputting  to  a teletype  or  line  printer. 

Data  suitable  for  graphics  display  is  controlled  by  the  parameter 
LOWESTLEVDSK.  The  values  (1  thru  7)  of  this  parameter  have  the  same 
meaning  as  defined  above.  A value  greater  than  7 produces  only  the 
mean  matrices  and  summary  statistics  (unformatted)  on  all  the  matrices 
starting  from  integration  time  of  10  seconds  up  through  the  maximum  of 
640  seconds. 

In  line  2 of  Fig.  2,  ^ defines  the  type  of  multiplication  used  by 
the  correlation  algorithm.  If  only  polarity  of  the  signal  is  used  to 
define  the  input  signal  then  LC  is  set  to  2.  For  linear  processing, 

LC:=1  indicates  linear  processing  utilizing  all  eight  data  bits  in  the 
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samples.  If  LC:=3  then  both  processes  are  performed 


In  line  3,  the  value  of  K is  used  to  reduce  the  redundancy  in  the 
sampled  input  data.  Only  every  Kth  sample  of  the  input  data  stream  is 
used  for  computation  of  the  covariance  matrices  and  statistics. 

In  line  4 and  5,  MTAU  (time  register)  and  QDLT  (doppler)  define  the 
resolution  of  the  covariance  computations  and  matrices.  The  values  are 
expressed  as  absolute  values.  In  terms  of  a matrix  this  refers  to  the 
rows  and  columns  on  either  side  of  the  center  lines.  + 16  MTAU  and 
+ 16  QDLT  thus  defines  a 33  x 33  matrix.  These  are  typical  values. 

In  line  6,  the  parameter  NTl  is  the  number  of  samples  used  by  the 
program  for  the  lowest  level  calculation.  It  is  equal  to  1/K  times  the 
number  of  samples  in  a level  1 (10  seconds)  integration  time,  (i.e., 

1/K  X (33  X 8)  samples/sec  x 10  secs).  For  the  program  in  Fig.  2 this 
gives  a value  of  83  for  NTl. 

In  line  7,  NTDNTl  defines  the  limits  of  integration  time  for  a 
computation.  It  is  the  ratio  of  the  number  of  samples  used  for  the 
longest  integration  time  to  the  smallest.  (Typically  640  seconds/10 
secs  * 64). 

In  line  8,  NTAU  is  the  incremental  shift  (in  units  of  data  samples) 
in  the  time  domain  for  each  computation  of  the  covariance  coefficient. 

In  line  9,  DELDEL  is  the  fundamental  doppler  increment  used  in  the 
covariance  coefficient  calculation.  It  is  a decimal  number  expressed 
as  a fraction  and  an  exponent  of  base  10. 

In  line  10,  is  a bias  value  associated  with  the  input  data.  If 
it  is  not  desired  to  start  the  calculations  until  a specific  time,  this 
value  Identifies  a starting  point  referenced  in  samples  from  the 
beginning  of  the  data. 

Lines  11  through  40  consist  of  15  pairs  of  parameters  for  each 
combination  of  processing  for  the  5 data  streams.  If  it  is  desired  to 
process  only  certain  pairs  of  data  then  only  those  parameters  associated 
with  that  pair  are  Included.  The  integers  in  brackets  indicate  which 
pair  of  data  streams  are  being  processed.  For  example,  in  lines  13  and 
14,  it  is  desired  to  crosscorrelate  received  data  set  2 with  the  set  1. 

A value  of  1 indicates  the  data  was  the  original  signal  used  as  a trans- 
mitting source  in  the  experiment.  The  values  2 through  5 are  receiving 
sites.  The  value  of  NYXA  for  a given  X,Y  pair  is  the  anticipated  time 
arrival  difference  (measured  in  samples)  between  the  two  recording  sites. 
The  value  of  DYXA  is  the  anticipated  doppler  shift  between  the  arriving 
signals. 
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Control  of  Experimental  Parameters 


Execution  of  a program  on  the  Illiac  IV  is  presently  done  by  batch 
processing.  Control  statements  are  stored  upon  the  disk  storage 
mechanism  of  the  Tenex  system  and  placed  in  a queue  for  later  execution 
by  the  Illiac  IV  complex.  A list  of  these  control  statements  is  shown 
in  Fig.  3.  When  these  statements  are  entered  into  the  batch  queue,  they 
become  known  as  the  Primary  Input  File  (PIF).  The  PIF  establishes  the 
proper  control  of  the  Illiac  IV  subsystems  and  directs  the  flow  of  data 
at  the  required  times. 

When  the  PIF  reaches  the  top  of  the  queue,  the  file  containing  the 
PIF  commands  is  opened,  and  the  control  statements  (known  as  ACL  state- 
ments) are  processed  in  the  sequence  shown  in  Fig.  3. 

The  PIF  consists  of  six  separate  functions  that  (1)  acquire  new 
parameters  for  the  correlation  algorithm,  (2)  combine  various  sub-modules 
of  the  SCI  correlation  program,  (3)  format  the  Illiac  disk  for  acceptance 
of  the  input  data,  (4)  move  data  from  magnetic  tape  storage  to  the  Illiac 
disk  (not  the  same  as  the  Tenex  disk  where  the  PIF  is  stored),  (5)  execute 
the  Illiac  processor  program  and,  (6)  place  the  program  results  on  the 
Tenex  disk  memory  where  it  can  be  transmitted  to  NRL. 

(1)  Input  New  Parameters 


On  each  request  for  a new  data  run,  various  parameters  can  be  mani- 
pulated in  the  correlation  algorithm  to  produce  the  desired  statistics 
and  covariance  coefficients.  As  the  program  is  presently  configured, 
these  parameters  are  placed  in  a special  program  that  is  compiled  by  the 
B6700  computer  and  then  placed  in  a Tenex  file  for  inclusion  into  the 
Illiac  when  the  NRL  request  for  a program  run  reaches  the  top  of  the 
queue.  At  the  present  time,  the  queue  is  executed  on  a first  in,  first 
out  basis  during  the  evening  hours  at  the  Illiac  complex.  If  a request 
by  NRL  is  placed  in  the  queue  during  the  day,  the  run  will  be  executed 
the  same  evening  and  be  availahle  for  inspection  and  display  by  the 
following  morning. 

The  parameters  that  may  be  manipulated  were  described  in  another 
section  and  shown  in  Fig.  2.  Once  these  parameters  are  corrected  as 
desired,  they  are  placed  in  the  Tenex  file  PARMVIWLA. SOURCE.  In  the 
overall  sequence  of  processing,  this  corresponds  to  step  lA  in  Fig.  1. 

Most  of  the  modules  used  by  the  correlation  function  are  written  in 
GLYPNIR.  This  is  a high  level  language  similar  to  FORTRAN  and  ALGOL. 
Special  features  are  included  that  are  unique  to  the  Illiac  processor. 

The  parallelism  required  to  perform  an  algorithm  efficiently  on  the  Illiac 
is  not  solved  by  GLYPNIR.  This  procedure  must  be  specified  in  detail  by 


the  programmer.  (In  Norse  mythology,  the  gods  despaired  because  a 
monster  roamed  uncontrolled  throughout  the  land.  Dwarfs  were  commissioned 
to  smith  a magic  chain,  called  GLEIPNER,  to  bind  the  monster  to  his  lair.) 

The  PIF  of  Fig.  3 takes  the  parameter  file  PARMVlWLA. SOURCE  and 
compiles  it.  The  control  command  in  the  PIF  that  performs  this  task  is 
shown  in  Fig.  3,  line  1.  The  command  GLYP  in  the  PIF  invokes  the 
GLYFNIR  compiler  which  resides  in  the  Borroghs  B6700  computer. 

PARMVlWLA. SOURCE  tells  the  compiler  the  name  of  the  Tenex  file  (Fig.  2) 
vrtiere  Che  newly  created  parameters  may  be  found.  With  this  information, 
the  compiler  proceeds  to  produce  another  file,  based  on  the  source  file, 
that  is  written  in  machine  language  acceptable  to  the  Illiac  processor. 
This  new  file  is  placed  on  the  Tenex  disk,  and  given  the  name  associated 
with  the  second  parameter  in  line  1 of  Fig.  3.  , PARMVlJF. REL  .(Parameters 
are  delimited  by  commas.) 

The  third  parameter  in  line  1,  PARMVlJF. LIST,  is  the  name  of  yet 
another  Tenex  disk  file.  In  this  file  a listing  of  compiled  statements 
in  Fig.  3 are  placed.  Any  mistakes  detected  in  the  source  file  by  the 
compiler  are  denoted  in  this  list  file. 

In  addition  to  GLYPNIR,  there  is  another  language  available  for 
programming  the  Illiac.  This  is  known  as  the  ASK  assembler.  ASK  is  a 
more  basic  language  specially  suited  for  programming  the  individual 
Processing  Elements  (PE's)  of  the  Illiac.  ASK  includes  the  instructions 
such  as  addition,  multiplication,  testing,  etc.,  that  are  part  of  the 
PE's  repertoire.  In  line  10  of  Fig.  3,  the  module  ASKVlJF  is  written 
in  ASK  assembly  language. 

(2)  Link  Modules  for  Execution 


The  compiler  output  files  (such  as  PARMVlJF. REL  described  above)  are 
formatted  in  a special  binary  machine  language  that  allows  the  program 
to  be  "relocated"  to  any  location  in  the  computer  memory.  Thus  these 
files  are  denoted  by  the  extension  REL  on  the  file  name.  It  is  the 
purpose  of  LINKED  (a  program  residing  in  the  Tenex)  to  take  all  the 
submodules,  "link"  them  together  as  one  complete  program,  and  produce 
a binary  file  which  is  no  longer  relocatable  but  assigned  to  absolute 
locations  in  Illiac  memory.  This  is  accomplished  in  lines  2 thru  11  of 
Fig.  3. 


The  acronym  ISV,  line  3,  instructs  LINKED  to  "save",  for  later 
execution,  all  the  files  in  the  absolute  program  that  will  later  be 
executed  on  the  Illiac  array.  Tlie  files  to  be  "included"  are  shown  in 
lines  4 thru  10.  These  included  files  are  documented  in  Ref.  1. 

Each  module  included  performs  a separate  function  on  the  Illiac  system. 
Note  that  in  line  8,  the  parameter  file  produced  in  step  1 is  included 
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in  Che  complete  program  to  be  linked  into  the  absolute  file  and  saved 
on  CASEVl.ISV, 

Line  11  of  the  FIF,  terminates  the  LINKED  program. 

(3)  Mapping  the  Illiac  Disk 

Before  any  data  transfers  can  actually  begin  on  I4DM,  a special 
program  called  MAP  must  be  invoked  to  properly  format  the  areas  on  the 
disk  memory.  The  MAP  subsystem  is  called  from  the  PIF  in  line  12.  In 
line  12,  MAPOUT  is  the  name  of  a Tenex  file  that  will  contain  allocation 
tables  used  to  establish  physical  areas  on  the  multidisk,  multihead 
I4DM  memory  system. 

The  I4DM  can  be  physically  visualized  as  containing  52  identical 
bands  (13  synchronized  disks  with  4 bands  on  each).  Each  band  is 
divided  into  300  contiguous  pages,  with  each  page  containing  1024  words 
of  64  bits  per  word.  The  page  is  the  smallest  area  that  can  be  accessed 
(read  or  write)  at  a given  transfer.  It  is  equivalent  in  length  to  one 
PE  (Processing  Element)  of  the  Illiac  IV  processor. 

The  allocation  tables  generated  by  MAP  are  created  from  FORMAT 
statements,  lines  13-22  of  Fig.  3,  that  define  data  areas  from  a pro- 
gramming standpoint.  This  greatly  simplifies  programming  efforts  since 
the  physical  considerations  mentioned  above  need  not  be  considered.  The 
allocation  tables  produce  a collection  of  pages  in  a specified  relation- 
ship to  each  other. 

The  FORMAT  control  statement  is  used  to  specify  an  area  on  the  I4DM 
and  associate  a name  to  it  that  may  be  referenced  by  the  program  within 
the  Illiac  processor.  In  lines  13-22  of  Fig.  3,  these  areas  have  been 
given  the  names  FZ,  FA,  FB,  FC,  FD,  FE,  and  FF.  The  quantity  shown  in 
parenthesis  after  each  area  is  called  a FORMATSPEC  that  identifies  the 
logical  location  and  length  of  each  named  area. 

There  are  two  main  commands  used  in  the  FORMATSPEC,  These  are  the 
skip  command  (denoted  by  S)  and  the  page  command  (denoted  by  P) . By 
use  of  these  operators,  the  entire  I4DM  may  be  formatted  according  to 
the  needs  of  the  program  and  the  length  of  the  input  and  output  data. 

The  MAP  subsystem  contains  a logical  pointer  that  may  be  positioned 
by  pages  anywhere  on  the  I4DM.  By  use  of  the  S operator,  the  pointer 
can  be  moved  forward  or  backward  n pages  from  its  present  position  by 
the  command  + n S.  At  the  location  of  the  pointer,  n pages  can  be 
reserved  by  the  command  + n P.  (A  negative  sign  would  reserve  pages 
from  the  pointer  backward  against  the  rotation  of  the  disk).  The  P 
command  resets  the  pointer  according  to  the  number  of  pages  reserved. 


8 


In  lines  13-19  of  Fig.  3,  various  page  lengths  and  areas  have  been 
reserved  for  the  five  data  files  that  will  contain  the  magnetic  tape 
data  from  the  four  receiving  sites  and  the  source  signal  recorded 
during  the  experiment. 

The  source  data  will  be  placed  in  the  area  FF.  From  line  15,  we  see 

that  this  area  starts  at  page  location  7201  on  the  I4DM  and  consists  of 

94  pages.  In  line  16,  the  pointer  is  positioned  7295  pages  into  the 
disk  area  with  90  pages  reserved  for  area  FB.  Note  that  area  FB  starts 
on  the  next  available  page  after  area  FF  (7201  + 94  = 7295).  It  is  not 
necessary  that  areas  be  contiguous  in  this  fashion,  but  it  is  good 
practice  to  minimize  the  total  space  used  on  the  I4DM. 

In  line  14,  one  page  at  the  beginning  of  the  source  data  is  reserved 

for  a label  identifying  the  particular  data  tape  and  run  number  that  is 
presently  being  processed.  At  first  it  might  appear  that  area  FA  is  95 
pages  long.  However  area  FF,  the  source  signal  data,  overlaps  area  FA 
except  for  the  page  location  7200.  Line  13  establishes  a relative 
starting  point  at  page  location  7200.  The  area  FZ  is  zero  pages  in 
length  and  is  used  to  establish  the  beginning  pointer  location  for  the 
Illiac  program. 

In  lines  20-22  space  on  the  I4DM  is  reserved  for  the  computed  results 
of  the  Illiac  program.  The  correlation  program  has  three  general  areas 
of  output  and  three  I4DM  disk  areas  are  reserved  for  each. 

In  line  22,  150  pages  have  been  reserved  for  the  packed  binary 
information.  This  includes  the  correlation  matrices  and  the  statistics 
(standard  deviation,  mean,  measure  of  significance,  peak  value)  that  are 
associated  with  each  matrix.  It  is  this  file  that  is  used  by  NRL  for 
producing  the  graphical  interpretation  and  printing  of  statistics  shown 
in  Figs.  5B,  6,  7,  8,  9.  The  actual  packing  of  these  values  is  described 
more  fully  in  another  section  and  shown  in  Fig.  4.  A complete  description 
may  be  found  in  Ref.  3. 

The  area  FOUT  contains  debugging  Information  produced  during  the 
execution  of  the  program.  In  the  event  of  failure  of  execution,  this  area 
may  be  examined  to  locate  any  problems  encountered  by  the  executing 
program. 

The  file  PRINTF  contains  information  similar  to  AREAO.  It  is 
formatted  in  ASCII  notation  so  that  it  may  be  typed  directly  on  a 
terminal  and  the  computed  values  may  be  examined  directly. 

Line  23  ends  the  MAP  subsystem.  The  configuration  tables  necessary 
to  completely  format  the  I4DM  are  now  complete  and  reside  on  the  file 
MAPOUT.  It  should  be  noted  that  the  MAP  subsystem  has  not  actually 
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required  the  I4DM  during  this  process.  Other  conanands  in  the  MAP 
subsystem  could  be  used  at  this  time  to  display  how  the  I4DM  would  look 
before  actual  program  execution  is  initiated.  In  this  way  errors  could 
be  detected  and  corrected  before  actual  execution  on  the  Illiac  array 
itself.  More  information  on  these  commands  as  well  as  the  entire 
Illiac  IV  control  language  can  be  found  in  Reference  4. 

In  line  24,  the  physical  disk  areas  of  I4DM  are  actually  reserved 
and  allocated  according  to  the  procedure  described  above.  The  ALLOC 
command  requests  use  of  the  I4DM  and  uses  the  first  parameter  MAPOUT 
(this  allocation  table  was  established  in  line  12  and  filed  with  inform- 
ation by  lines  13-22)  as  an  identifier  to  find  the  allocation  tables 
on  the  Tenex  disk. 

The  second  parameter,  ALGID,  is  a name  that  will  be  assigned  to  all 
the  areas  that  are  allocated  to  this  particular  job.  During  the  period 
that  this  identifier  is  in  effect,  no  other  users  of  the  Illiac  processor 
can  occupy  the  spaces  reserved  on  the  I4DM  under  this  name.  When  the 
job  is  completed,  this  name  will  be  used  to  deallocate  the  reserved 
areas  and  allow  their  use  by  other  jobs. 

(4)  Loading  the  IV  Disk 


We  are  now  in  a position  to  begin  data  movement  onto  the  I4DM  disk. 
This  is  accomplished  with  the  MOVE  instructions  in  lines  25-30.  The 
procedure  is  to  acquire  data  from  the  Tenex  disk  and  move  it  to  the 
14dm  disk. 

Since  there  were  five  sites  recorded  during  the  field  experiment, 
five  files  of  data  are  transferred  to  the  14dm  disk.  This  data  was 
obtained  from  magnetic  tape  and  placed  on  the  Tenex  disk  before  the  PIF 
was  entered  into  the  Illiac  queue.  The  five  sites  are  identified  in 
the  first  parameter  of  the  MOVE  commands.  In  lines  26-30,  these  data 
files  are  Identified  as  run  41C2,  sites  SO  thru  S4. 

These  data  files  are  assigned  to  the  I4DM  in  the  second  parameter  of 
the  MOVE  command.  Note  that  these  names  (FE  thru  FF  in  lines  26-30)  are 
identical  to  the  areas  that  were  allocated  by  the  MAP  subsystem  in 
lines  15-19. 

Line  25  is  a label  file  used  for  identification  of  the  magentic  tape 
data.  As  explained  before,  it  consists  of  1 page  of  information.  The 
programs  executed  at  NRL  access  this  information  to  identify  the 
particular  data  run. 

(5)  Executing  the  Illiac  Program 

Line  31  is  the  line  we  have  all  been  waiting  for.  Until  now  we  have 
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only  been  making  preparations  to  execute  a program.  The  new  parameters 
have  been  compiled  in  line  1 of  the  PIF.  The  programs  have  been  linked 
together  in  lines  2 thru  11;  and  the  Illiac  IV  disk  has  been  allocated 
and  data  entered  into  the  proper  areas  in  lines  12  thru  30. 

The  RUN  command  enters  the  programs  (linked  and  identified  as 
CASEVl.ISV  in  line  3 of  the  PIF)  into  the  Control  Unit  (CU)  and  the 
Processing  Elements  (PE's)  of  the  Illiac  IV  processor.  As  the  program 
is  executing,  it  will  request  input  data  from  the  I4DM  by  using  the  area 
names  FA,  FE,  FD,  FC,  FB,  FF.  As  the  correlation  coefficients  are 
computed,  they  are  written  to  the  I4DM  using  the  areas  PRINTF,  FOUT, 
AREAO. 

The  second  parameter  of  line  31  allocates  the  maximum  time,  in 
seconds,  that  the  Illiac  IV  processor  can  use  for  this  particular  run. 

In  the  unlikely  event  that  the  program  loops  continuously,  this  parameter 
will  terminate  the  run  at  the  allocated  time. 

(6)  Saving  the  Computed  Results 

When  the  program  has  finished  executing,  the  batch  queue  processor 
goes  to  the  next  line,  32,  to  identify  the  areas  on  the  Tenex  disk  that 
will  be  used  for  storage  of  the  computed  results.  As  explained  before, 
these  areas  are  PRINTF,  FOUT,  AREAO.  In  the  particular  PIF  shown  in 
Fig.  3,  only  the  binary  data  file,  AREAO  is  permanently  saved  on  the 
disk.  This  is  the  file  used  by  NRL  for  formation  of  the  graphic  data 
and  compilation  of  statistics  described  in  another  section  and  shown  in 
Figs.  5B,  6,  7,  8,  9.  As  shown  in  line  32  of  the  PIF  this  saved  area 
is  given  the  Tenex  file  name  0UTV2.DISK.  If  the  areas  FOUT,  PRINTF  are 
to  be  saved,  additional  MOVE  statements  would  be  Included  in  the  PIF 
with  Tenex  file  names  to  identify  them. 

Line  33  terminates  the  PIF  by  deallocating  the  I4DM  reserved  areas. 
These  areas  are  now  free  for  the  next  user  in  the  batch  queue.  The 
next  sections  describe  how  the  information  from  the  Illiac  IV  data  is 
extracted  for  transmission  and  processing  at  NRL. 

NRL  Data  Acquisition  and  Display 

The  binary  file  displayed  in  Fig.  4 contains  tabulated  statistics 
and  the  covariance  coefficients  computed  for  each  matrix.  NRL  has 
developed  algorithms  for  unpacking  the  data  from  the  Tenex  disk  file 
and  transmitting  it  to  NRL.  Additional  algorithms  transform  the  correla- 
tion coefficients  to  display  vectors  suitable  for  graphic  display  as 
a simulated  3-D  plot  with  isometric  projection  and  hidden  line  removal. 
For  examination  of  actual  parameters,  the  statistics  may  be  collated 
according  to  a particular  site  pair  and  tabulated  for  each  matrix  in 
tht  computation. 
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Graphics  Display 


Graphics  display  of  the  data  requires:  (1)  accessing  the  binary 

data  file  containing  the  computed  covariance  coefficients,  (2)  unpacking 
a matrix  from  the  binary  file  via  parameters  Issued  from  the  NRL 
terminal,  (3)  transmission  of  the  matrix  coefficients  to  NRL,  (4)  forma- 
tion of  a simulated  3-D  display,  (5)  Identification,  extraction  and 
display  of  peak  Information.  As  shown  In  Fig.  4,  all  matrices,  along 
with  Identifying  statistics,  are  packed  within  one  Tenex  disk  file.  In 
order  for  the  unpacking  algorithm  to  access  a desired  matrix,  certain 
minimum  information  is  required  by  the  program.  The  program  is  written 
to  query  the  operator  interactively  as  it  requires  information  to  select 
a particular  matrix.  In  Fig.  5A,  the  program  has  been  loaded  and 
immediately  requests  the  NRL  operator  to  identify  the  particular  Tenex 
file  that  is  to  be  displayed.  In  Fig.  5A,  the  underlined  alphanumerics 
a.re  the  questions  the  program  asks  of  the  operator.  The  operator 
responds  (not  underlined)  until  sufficient  information  is  supplied  to 
identify  a unique  matrix. 

For  the  example  shown,  the  operator  has  requested  the  file 
0UTV2.DISK.  The  name  of  this  file  was  created  during  the  batch  run  on 
the  Illiac  IV.  This  name  was  allocated  by  line  32  of  the  Primary  Input 
File,  Fig.  3. 

In  the  next  line,  the  program  requests  the  operator  to  confirm 
that  output  data  is  to  be  sent  to  the  terminal  that  loaded  the  program. 
This  is  accomplished  by  typing  TTY:  and  confirming  with  a carriage 
return.  In  the  future  it  is  expected  that  more  than  one  terminal  or 
data  path  may  be  active  at  a time.  Typing  in  a different  acronym  hare 
would  divert  the  output  as  desired. 

The  program  then  requests  the  operator  to  identify  the  run  number 
that  is  presently  being  processed.  This  information  is  used  to  identify 
the  graphics  display.  These  run  numbers  correspond  to  the  notation  used 
during  the  experiment  to  identify  the  type  of  signal  and  its  origin. 

On  the  next  two  lines,  the  program  requests  the  operator  to  identify 
the  signals  that  were  correlated.  In  this  particular  example,  site  0 
(the  source  signal)  is  crosscorrelated  with  site  1. 

In  the  next  line,  the  operator  selects  which  type  of  processing  he 
desires.  In  the  matrix  shown.  Fig.  5B,  linear  was  selected.  Both 
Linear  (8  bits  of  magnitude)  and  Clipped  (sign  only)  versions  of  the 
data  can  be  made  available  in  the  same  binary  file.  This  parameter  was 
established  in  line  l of  Fig.  2 during  the  batch  execution  of  the  Illiac 
programs.  As  shown,  and  discussed  elsewhere,  the  value  3 instructs  the 
program  to  produce  clipped  and  linear  versions  of  the  covariance  matrices. 
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Finally,  the  operator  must  establish  the  level  of  integration  time 
and  the  interval  within  that  integration  time.  As  shown  in  Fig.  5 A,  an 
integration  time  of  6 is  selected.  Integration  times  can  assume  the 
values  10  sec,  20  sec,  40  sec,  80  sec,  160  sec,  320  sec,  and  640  sec. 

The  LEVEL  1 requested  in  Fig.  5A,  selects  an  integration  time  of  10  secs. 

A level  7 request  would  yield  the  matrix  that  corresponds  to  an  inte- 
gration time  of  640  seconds. 

' The  value  of  RK  identifies  which  interval  in  a given  integration  time 

, is  requested.  Since  level  7 is  the  highest  integration  time,  there  is 

only  one  interval  at  level  7.  Therefore  a value  of  0 is  selected.  j 

(The  numbering  system  runs  from  0 to  63  and  not  1 to  64).  For  LEVEL  1, 
there  are  64  ten  second  intervals  available.  Should  an  invalid  value  be 
t typed,  the  program  will  type  an  error  message  and  request  a new  value. 

The  values  of  RK  and  LEVEL  are  more  easily  vistialized  in  tabulations 
such  as  Fig.  6,  where  each  matrix  for  a given  RK  (the  leftmost  column) 
is  tabulated  against  the  seven  integration  times,  Tl  thru  T7. 

At  this  point,  all  information  needed  to  identify  a unique  matrix 
has  been  introduced  and  the  graphics  program  exists  from  the  interactive 
mode  and  begins  unpacking  the  binary  data  file.  As  the  matrix  is 
unpacked,  the  header  (described  in  another  section  and  Fig.  4)  is 
examined  and  displayed  on  the  CRT,  along  with  the  graphic  vectors,  as  a 
confirmation  of  the  identity  of  the  matrix.  As  shown  in  Fig.  5B,  the 
header  indicates  this  matrix  to  be  the  results  of  site  pair  0,  1 for 
run  41C2.  It  is  a linear  matrix  with  LEVEL=6  (320  sec)  as  the  operator 
requested  in  Fig.  5A.  The  other  statistics  shown  in  Fig.  5B  further 
identify  the  matrix  and  the  parameters  that  produced  it.  These  parameters 
were  described  in  the  section  dealing  with  Fig.  2. 

At  the  top  of  Fig.  5B  are  two  cross-plots  that  precisely  display 
the  peak  value  of  the  plot.  The  display  program  searches  the  matrix 
for  the  largest  value  and  then  extracts  the  row  (TAU  or  time  dimension) 
and  the  column  (DELTA  or  doppler  dimension)  associated  with  the 
displayed  matrix. 

As  shown,  the  TAU  crossplot  displays  the  peak  in  the  hidden  line 
display  as  if  no  other  lines  were  present.  The  TAU  crossplot  shows  the 
precise  location  of  the  peak  in  relation  to  the  33  points  of  the  correla- 
tion function.  From  the  parameters  (Fig.  2)  introduced  into  the  Batch 
stream  when  the  Illlac  program  was  run,  the  exact  TAU  value  could  be 
ascertained.  (Automatic  labeling  of  the  tic  marks  with  the  actual 
values  of  TAU  and  DELTA  will  be  done  in  the  next  version  of  the  graphics 
program) . 

Similarly  a crossplot  along  a specific  column  is  produced  to  give 
better  insight  into  the  DELTA  parameter.  This  plot  looks  "into"  the 
display  at  the  column  number  containing  the  peak. 
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Statistical  Parameter  Tabulations 


For  the  correlation  of  a pair  of  recorded  signals,  the  statistics  of 
each  integration  interval  may  be  examined.  This  is  done  by  collating 
the  statistics  of  each  matrix  and  tabulating  them  individually.  This 
is  shown  in  Figs.  6,  7,  8,  and  9 where  the  means,  standard  deviation, 
measures  of  significance  and  peak  values  of  each  matrix  are  displayed. 

Once  the  program  is  executing,  information  is  requested  of  the 
operator  to  identify  the  statistics  desired.  The  procedure  involving 
the  operator  is  similar  in  nature  to  the  one  described  under  graphics 
display.  At  present,  the  statistics  program  described  here  and  the 
graphics  program  described  previously  are  separate  entities.  In  the 
near  future  these  programs  will  be  linked;  then  computer  algorithms, 
under  operator  control,  will  select  and  display  only  those  matrices 
that  are  of  interest  to  the  operator.  Additional  interactive  commands 
will  be  included  to  allow  the  operator  to  establish  constraints  upon  the 
algorithm  and  select  the  values  of  interest.  This  will  permit  the 
operator  to  selectively  extract  information  from  the  data  base. 

As  before,  the  requests  of  the  program  are  underlined,  with  the 
operator  response  following.  In  Fig.  6,  for  example,  the  Illiac  IV 
data  file  0UTV2.D1SK  is  identified  as  the  data  base  for  the  tabulations 
to  follow;  this  file  of  data  represents  Run  41C2  of  the  experiment  as 
shown  in  the  third  line  of  Fig.  6.  In  the  next  two  lines,  site  0 and 
site  1 are  the  site  pair  we  wish  to  examine.  This  is  followed  by  a 
declaration  to  select  only  linear  correlation  values.  Finally  the 
particular  statistic  desired  is  declared. 

Given  this  information,  the  program  searches  for  the  desired  statis- 
tics and  prints  out  a header  that  confirms,  from  information  contained 
in  the  binary  file,  the  input  parameters.  The  values  of  each  matrix 
are  printed  according  to  the  integration  time  and  the  interval  within 
the  640  second  data  base.  As  shown,  the  value  of  RK  along  the  left-hand 
side  determines  the  interval  number.  The  integration  times  are  T1  thru 
T7  representing  10  second  thru  640  second  integration  times.  As  can  be 
seen,  there  are  sixty-four  10  second  intervals,  thirty-two  20  second 
intervals,  etc.  At  level  T7  (640  seconds)  there  is  only  one  value  since 
this  correlation  integral  has  the  limits  of  the  entire  data  base. 

After  the  last  RK  values  are  listed,  the  mean  value  of  the  individual 
integration  times  is  listed. 

Figs.  7 and  8,  the  means  and  standard  deviations  of  the  matrices  are 
tabulated  in  a similar  manner  to  the  measures  of  significance. 

Fig.  9 shows  the  tabulations  for  the  peaks  and  the  indices  of  the 
peaks  for  each  matrix.  These  values  may  be  used  to  establish 
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calibration  values  on  the  graphic  display  of  Fig.  5B.  From  Fig.  5B,  we 
find  that  this  matrix  was  produced  from  Run  41C2,  sites  0 and  1,  with  an 
RK=0  interval  value  in  level  6.  From  Fig.  9,  we  find  that  this  particular 
matrix  (level  T6  and  RK=0)  has  a peak  value  of  5156.  A doppler,  or 
delta,  index  value  (DI)  of  18  is  associated  with  this  peak.  The  time, 
or  tau,  index  (TI)  has  the  value  16.  This  can  be  crosschecked  on 
Fig.  5B  from  the  crossplots  at  the  top.  The  tick  marks  are  numbered  1 
to  33  with  a value  of  17  in  the  center.  As  can  be  seen,  the  peak  occurs 
at  DELTA  and  TAU  values  18  and  16  respectively.  In  the  near  future,  the 
values  of  Fig.  9 will  be  used  to  automatically  label  the  diagrams  of 
Fig.  5B,  by  linking  the  two  programs  together. 

Illlac  IV  output  data 

The  Illiac  IV  program  assumes  that  Tenex  disk  memory  is  partitioned 
into  64  bit  words  with  8 bytes  per  word.  The  covariance  coefficients 
amplitudes  have  been  restricted  to  one  byte  in  length.  Thus  eight 
coefficients  can  be  neatly  packed  in  one  word.  A Tenex  page  is  defined 
as  2048  bytes. 

A covariance  coefficient  matrix  is  defined  as  S]^=(2  tau  + 1)(2  delta 
+ 1)  elements,  where  tau  and  delta  define  the  limits  of  the  matrix.  It 
is  anticipated  that  M and  Q will  remain  fixed  at  16  for  the  near  future. 
This  gives  a matrix  of  33  by  33  or  1089  distinct  elements.  As  the 
matrices  are  computed,  they  are  stored  sequentially  on  the  Tenex  disk 
memory.  In  order  to  identify  each  matrix,  a 96  byte  header  precedes  each 
matrix.  Thus  ST=(Si+96)*1185  bytes  are  required  for  each  matrix. 

By  convention,  each  matrix  and  header  start  on  a page  boundary.  If 
the  total  byte  count,  ST,  does  not  fill  the  complete  page  of  2048  bytes 
the  remainder  of  the  page  is  left  unused.  Fig.  4 shows  the  format  of 
the  header  and  the  matrix. 

Five  separate  signals  can  be  used  in  the  correlation  program.  As 
explained  in  another  section,  the  XY  pairs  that  are  used  as  input  data 
to  the  Illiac  correlation  program  can  be  selected  by  the  parameter  NYXA 
in  Fig.  2. 

The  XY  pairs  are  computed  and  stored  in  the  following  sequence: 


(0,0), (1,1), (2, 2). (3, 3), (4, 4), (0,1), (0,2), (0,3), (0,4), (1.2), (2, 3), 
(3, 4), (1,3), (2,4), (1,4). 


If  a particular  XY  pair  is  not  requested,  no  space  is  allocated  to 
it  in  the  file.  The  integers  0 thru  4 are  used  to  identify  the  sites 
associated  with  the  raw  field  data.  The  ship  source  signal  is  denoted 
by  0. 
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NRL  Terminal  System 


The  NR'f,  Terminal  System  is  shown  in  Fig.  10.  The  system  is 
configured  for  flexibility  in  acquisition  and  display  of  information 
available  on  the  ARPANET.  The  total  system  consists  of  a PDP-11/40 
central  computer,  an  Imlac  PDS-4  refresh  CRT,  a Versatec  llOOA  Printer/ 

Plotter,  and  a Computek  415  storage  screen  CRT. 

For  acquisition  of  ARPANET  data,  the  11/40  interfaces  to  a Bell 
System  202R  modem  with  a direct  line  to  a Terminal  Interface  Processor 
(TIP)  located  at  the  Mitre  Corporation,  McLean,  Virginia.  The  202R 
modem  allows  simultaneous  (full-duplex)  transmission  at  a 1200  bit/sec 
rate  between  the  NRL  system  and  the  ARPANET.  The  202R  uses  asynchronous 
transmission  (start  and  stop  bits  on  each  8 bit  data  character). 

As  shown  in  Fig.  10,  ARPANET  Information  is  received  and  sent  by  a 
full  duplex  module  in  the  PDP-11.  Also  shown,  in  dotted  lines,  is  a 

similar  connection  to  a full  duplex  module  in  the  Imlac  PDS-4.  With  i 

this  direct  connection  the  Imlac  CRT  could  receive  graphics  information 

directly  from  the  network  and  bypass  the  PDP-11  computer.  This  mode  i 

would  be  useful  for  direct  viewing  of  graphics  information  available  at  i 

the  host  computer.  This  direct  mode  of  operation,  however,  does  not  I 

allow  much  flexibility  in  operation  since  the  CPU  associated  with  the  • 

Imlac  is  designed  primarily  to  support  graphics.  This  would  place  a j 

much  larger  burden  on  the  host  computer  (the  I4-Tenex)  to  provide  the  j 

operator  interaction  and  data  manipulation  that  will  be  required  of  the 
NRL  system  to  effectively  analyze  the  large  data  base  produced  by  the 
Illiac  IV.  Responsive  operator  interaction  would  be  inhibited  in  this 
mode. 

$ 

The  output  of  the  Illiac  IV  consists  of  a data  structure  shown  in 
Fig.  4.  This  file  can  be  transmitted  in  its  entirety  to  the  PDP-11 
where  local  programs  can  direct  the  data  streams  along  several  paths 
under  operator  control.  In  most  cases,  the  entire  data  file  will  be 
transferred  to  magnetic  tape  for  permanent  storage  of  the  Illiac  IV 
computed  results.  For  on-line  operation  this  file  can  be  stored  on 
the  disk  memory  where  rapid,  random  access  to  the  individual  matrices 
and  statistics  are  available.  The  disk  is  used  in  all  operations  that 
support  graphics. 

The  matrix  shown  in  Fig.  5B  consists  of  1089=33x33  discrete  points. 

The  amplitudes  are  quantized  in  8 bit  samples.  If  each  sample  contains 
a start  and  stop  bit,  as  required  for  asynchronous  transmission,  then 
1089x10=10890  bits  must  be  transmitted  for  each  matrix.  At  a 1200  bit/ 
second  rate,  this  would  require  about  9 seconds  to  completely  transmit 
the  matrix  information  from  the  host,  I4-Tenex,  to  the  NRL  system.  In 
reality,  the  transmission  takes  longer  and  is  variable  dependent  on  the 

number  of  users  logged  onto  the  I4-Tenex  system,  the  amount  of  traffic  I 
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on  the  network  itself,  computer  overhead,  etc.  If  operation  near  this 
maximum  rate  is  desired,  the  data  transfer  can  be  done  during  the  less 
busy  hours  such  as  early  morning. 

The  present  disk  storage  capability  of  the  PDP-11  is  2.4  million 
bytes  (8  bits)  of  storage.  If  we  assume  half  of  the  disk  is  available 
for  data  storage  (the  other  half  is  reserved  for  the  monitor  and 
programs)  then  about  1100  matrices  (1.2  x 10®  bytes/1.089  x 10^  bytes/ 
matrix)  are  available  for  rapid  display  and  analysis  at  any  one  time. 

If  the  display  device  is  the  Imlac  PDS-4,  the  disk  data  can  be 
transmitted  at  a very  high  rate  of  transfer.  The  Imlac  CPU  connects 
directly  to  the  PDP-11  Unibus.  This  allows  data  transfers  at  memory 
speeds.  High  speed  transfers  are  essential  if  displays  require  rapid 
modification  (such  as  on-line  rotation)  and  continuity  of  the  CRT  image 
is  to  be  maintained  without  flicker.  Estimation  of  the  maximum  number 
of  images  (matrices)  that  can  be  transferred  between  the  two  devices 
is  difficult  since  applications  software  must  be  included  in  the  timing 
estimates.  Actual  data  transfers  would  be  in  the  multikilobit  range; 
overhead  to  support  data  acquisition  from  the  disk  and  programs  to  format 
raw  data  to  graphic  images  and  interact  with  operator  commands  will  add 
to  the  overhead  and  dilute  the  high  transfer  speeds.  For  realistic 
applications,  however,  this  will  pose  no  problem;  images  must  be  viewed 
by  operators, and  parameters  extracted,  for  each  image  displayed.  The 
actual  requests  for  new  data  by  the  operator  will  be  quite  slow  compared 
to  hardware  speeds. 

The  Imlac  PDS-4  contains  a unique  hardcopy  device  that  allows  any 
arbitrary  display  on  the  CRT  to  be  photographically  copied  to  sensitive 
paper.  Hardcopy  devices,  such  as  the  one  associated  with  the  Computek 
in  Fig.  10,  are  quite  commonly  used  for  storage  screen  CRT  devices. 

Since  the  image  is  statically  stored  on  a phosphor  screen  it  can  be 
easily  scanned  electrically  for  line-by-line  serial  transfer  to  sensitive 
paper.  Unlike  the  storage  screen  devices,  a refresh  CRT  always  displays 
a d}mamic  picture  from  vector  information  contained  in  memory.  This  | 

vector  information  cannot  be  scanned  serially  by  line-by-line  techniques. 

In  the  Imlac  hardcopy  device,  a separate  small  CRT  image  is  projected 
with  lenses  onto  sensitive  paper.  The  hardcopy  CRT  is  connected  in 
parallel  with  the  main  CRT;  thus  any  image  is  easily  copied. 

Under  its  present  configuration  the  Imlac  can  display  8 thousand 
vectors  of  arbitrary  length  within  the  screen  area.  If  the  vectors  are 
less  than  full  screen  size,  fewer  bits  are  needed  to  define  them  with 

the  result  that  even  more  vectors  can  be  included  in  a display.  j 

The  display  area  of  the  Imlac  is  divided  into  2048  x 2048  individually 
addressable  units.  There  are  three  vector  modes  for  addressing  a 
particular  point  on  the  screen:  (1)  Long  vector  mode  uses  two  words 
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(16  bit  words)  of  memory  to  define  a vector.  This  allows  maximum  "jumps" 
across  the  screen,  (2)  Medium  vector  mode  requires  one  word  of  memory. 

A maximum  of  64  display  units  can  be  spanned  in  this  mode,  (3)  In  the 
short  vector  mode  the  longest  movement  is  three  units.  These  vectors 
are  packed  two  in  a word. 

The  Imlac  PDS-4  also  contains  a scaling  instruction  that  multiplies 
the  above  vector  lengths  by  factors  up  thru  8.  This  allows  even  more 
flexibility  when  defining  vectors  and  permits  very  compact  storage  when 
proper  programming  techniques  are  applied. 

Additional  features  include  a hardware  "blink"  feature.  Any  figure 
displayed  can  be  made  to  blink  on  and  off.  The  Imlac  also  contains 
programmable  intensity  they  can  add  16  level  gray  tones  to  figures 
on  the  screen.  There  is  also  a hardware  rotation  figure  that  permits 
rotation  by  multiples  of  90°. 

Other  graphics  capability  in  Fig.  10  include  a Versatec  llOOA 
printer/plotter.  The  Versatec  uses  an  electrostatic  writing  principle 
that  uses  line-by-line  writing  techniques.  Each  line  consists  of 
writing  nibs  packed  on  .01  inch  centers.  Characters  or  graphics  are 
produced  by  breaking  the  images  into  separate  full  length  (about  10 
inches  on  the  llOOA)  scans  with  each  nib  containing  an  electrostatic 
charge  proportional  to  the  "slice"  of  the  image  that  is  being  displayed. 
The  nibs  electrify  the  moving  paper  which  passes  over  a liquid  toner 
containing  suspended  carbon  particles  which  adhere  to  the  charges  on  the 
paper.  In  plotting  mode,  there  are  120  scans  per  second  with  the  paper 
supply  moving  at  1.2  inches  per  second.  Thus  a 10  inch  x 10  inch  plot 
can  be  completed  in  12  seconds  with  a 100  point/inch  resolution  in  both 
the  X 6e  Y directions. 

Additional  features  allow  the  plotting  mode  and  character  printing 
mode  to  operate  simultaneously,  permitting  alphanumeric  annotation  of 
any  graphics  display.  In  the  alphanumerics  mode  there  are  132  columns 
of  characters  per  line  and  6.6  printed  lines  per  inch.  Maximum  printing 
speed  is  500  lines  per  minute. 

The  NRL  system  includes  both  seven  and  nine  track  magnetic  tape 
facilities.  Nine  track  tape  units  are  used  primarily  for  the  experimental 
data  since  the  computed  correlation  coefficients  are  quantized  to  eight 
bits  of  magnitude.  The  seven  track  unit  is  available  for  compatibility 
with  the  source  signal  recordings  of  the  experiment.  Additionally,  this 
seven  track  unit  provides  access  to  previously  recorded  data  bases 
collected  by  NRL  during  field  operations  of  other  experiments. 
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Fig.  2 — Illiac  IV  parameter  list 
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I GLY.->  JF.  Lli  r 

9  LIMK£D 

3 isy  CASEv;i . iSsJ 

4 INCLJDE  MAINv/1  JF.KEL 

5 INCLJDE  DliKJl  JF.HEL 

6 INCLJDE  G INI  rjlJF..iEL 

7 INCLJDE  INf^jrjD./iEL 

8 INCLJDE  PAAMJIJF.HEL 

9 INCLJDE  DELTAJl JF.HEL 

10  INCLJDE  AJKJIJF.HEL 

I I END 

12  MAF  ,MAF>0Ur 

13  FORMAT  FZ,  (72003  #0r>  ) 
lA  FOfiMAT  FA,  ( 72003  >95F) 

15  FOilMAr  FF, (72013, 94j3) 

16  FOilMAI  FB, (72953, 90P) 

17  FO.-iMAr  FC,(73  853,90F>) 

18  FOHMAT  FD,  (74  753,90?) 

19  FOBNAT  FE, (75653,90?) 

20  FOrfw^^f  ?.-iINTF,  (77003,100?) 

21  FOtiMAT  FOOT,  (78003, 10?) 

22  FOiiNAT  ArlEAO,(  79003, 150?) 

23  END 

24  ALLOC  MA?OJT,ALCID 

25  MOJE  LAaL£-33.41G2, 14DM:FA 

26  MOVE  fiJN-33.4lC2-30,  I4DN:FE 

27  N07E  tiJN-33.41C2-31  ,I4DN:FD 

28  MOOE  rlJN-3  3.41C2-32,  I40N:FC 

29  MOJS  HJN-33.41C2-33,I4DM:F3 

30  NOJS  3GI-17.41C2-34, 14DM:FF 

31  AUN  CA3Syi . 13 J,NAXTIN=3600 

32  NOJE  I4DM ; AHEAD,  rJT”G. 2 1 E” 

33  DALLOG  ALGID 

Fig.  3 — Primary  input  file  (PIF) 
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[BEGINNING  OF  BINARY  FILE- 
OFFSETS  DENOTED  IN  BYTES  18  BITS) 
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Fig.  5(a)  — Operator  matri.x  selection 


< .vL  H > ) J 1 -J  i • J I j L \L  A 


fit,-:  VAV£  o.-r  ILLIVJ  J\i\; 

V 1 •*:  I r I ; lit:  [.j  I )\i » I .ir-i  j 
I' iv-:  tjv  MjvB--:.;;  j 1 :,2 


jlf£  X=  0 
all'-:  t=  1 

OMi^r>-;i3  0.1  LI  V-:  \.i 

i'i'r>-:  ^XAin.t'  of  rn-:  --ollovimiI 
■*s 3 V i il'007  r*-^  \'Ko  I i J I 1' 

■^s  iTm 


a 1'  4 1 I 

a 1 UAL 

r’A.t  ACHll  la 

■t  1 At3  JL  A 1'  I Oi\  j 

4U2 

3 I l ^ A = . 

) 

ill'-:  f=  1 

LIM£A.t 

3 3 r(4. 

u 

0?  i 

I IMIF I J AMO A 

■ tK 

f 1 

fi 

1.3 

i'/4 

1'  6 i'  6 

1) 

/I  . .4 

4.6 

4.7 

6.9  10.7 

1 

4 .A 

4 . 4 

5.3 

5 • S 

7 . d 10.6 

j 

4.6 

4.4 

4.6 

'5.  J 

7. 1 

) 

4.4 

4.7 

4.7 

7.0 

7.d 

/4 

4 . 4 

4 . 6 

6.0 

3.7 

5 

4.4 

4.7 

4 . d 

o • 0 

<3 

4 . • 3 

4.6 

4 . 9 

6 • *:i 

7 

4.4 

4.6 

6 . 1 

6*0 

6 

4 . 5 

4.6 

6 . 1 

9 

4.7 

4.6 

4 . 7 

1 ) 

4 . 6 

4 . d 

4. 

1 1 

4.7 

4.6 

4.  ; 

1 J 

4 . S 

4.6 

6.1 

1 I 

4.4 

4 .4 

6.1 

4.4 

4.« 

4.4 

l-> 

4 . 4 

4 . d 

4.6 

16 

4.4 

4.6 

17 

4.  1 

4.4 

16 

4.6 

4.3 

l ) 

4.6 

4.6 

n 

4.7 

4.3 

n 

4 .d 

4.d 

Fig.  6 — Measures  of  significance  tabulations  (Continues) 
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Fig.  10  — NRL  compulor  systom 


